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The stability of the nuclear matter system with respect to density fluctuations is examined 
exploring in detail the pole structure of the electro-nuclear response functions. Making extensive 
use of the method of dispersion integrals we calculate the full polarization propagator not only 
for real energies in the spacelike and timelike regime but also in the whole complex energy plane. 
The latter proved to be necessary in order to identify unphysical causality violating poles which 
are the consequence of a neglection of vacuum polarization. On the contrary it is shown that 
Dirac sea effects stabilize the nuclear matter system shifting the unphysical pole from the upper 
energy plane back to the real axis. The exchange of strength between these real timelike collective 
excitations and the spacelike energy regime is shown to lead to a reduction of the quasielastic peak 
as it is seen in electron scattering experiments. Neglecting vacuum polarization one also obtains 
a reduction of the quasielastic peak but in this case the strength is partly shifted to the causality 
violating pole mentioned above which consequently cannot be considered as a physical reliable result. 
Our investigation of the response function in the energy region above the threshold of nucleon 
anti-nucleon production leads to another remarkable result. Treating the nucleons as point-like 

1.^ ' Dirac particles we show that for any isospin independent NN-interaction RPA-correlations provide 

J^ . a reduction of the production amplitude for pp-pairs by a factor 2. 



> 

Cn , I. INTRODUCTION 

o 

*vj • Dirac sea effects are perhaps the most interesting characteristic of relativistic many body theory. In nuclear 

f~^ ■ matter e.g. the effective mass of the nucleons and their coupling to the mesons undergo considerable corrections 
^\ , due to vacuum fluctuations M, 0. In addition, relativistic many body theory predicts strong corrections of the 
0^ properties of hadrons at high nuclear matter densities. Unfortunately, up to now it proved to be very difficult to 
^^ , find unambiguous experimental evidence for such effects. One method to probe the Dirac sea is to investigate the 
-*— > ■ production of antiprotons and other antibaryons in heavy ion collisions. Such experiments are currently performed by 
several groups under quite different kinematic conditions. These range from subthreshold production of antiprotons 
in Ni+Ni collisions at 1.85 GeV/nucleon at SIS (GSI) up to antibaryon production in Au+Au collisions at 14.6 
GeV/nucleon at AGS (E878 Collaboration at BNL) g] and sulphur-nucleus collisions at 200 GeV/nucleon at SPS 
(NA35 collaboration at Cern) . Beside the search for evidence of the quark-gluon plasmaj || the attention is drawn 
to possible changes of the properties of the nucleons and especially the vector mesons R] - p| in the nuclear medium. 
The presence of a filled Dirac sea of antiparticles, in addition to the Fermi sea of protons and neutrons, has 
^ ■ consequences not only for the static properties mentioned above, but also for the dynamics of nuclear matter. The 
C^ \ latter may be probed e.g. by measuring the response of the system to the disturbance induced by an external virtual 
photon, as it is done in inelastic electron scattering experiments. A well known and as yet unsolved problem in this 
field is an apparent quenching (in comparison with the prediction of simple one-particle models) of the longitudinal 
contribution to the inclusive scattering cross section, particularly for medium heavy and heavier nuclei like ^''Ca, 
^^Fe and ^^^U. First observed for momentum transfers of some lOOMeV/c ||l^, Q, a SLAG experiment [^ later 
showed that the quenching even persists up to 1 GeV/c. This may be taken as a hint that relativistic effects play 
a role. Further, it seems reasonable that in heavier nuclei collective phenomena reduce the strength of longitudinal 
one-particle excitations. So the quenching problem of the longitudinal response function appears as an appropriate 
subject for the application of relativistic many body theory. 

The first conserving approximation which also includes collective phenomena in a proper way is the random phase 
approximation (RPA). RPA-calculations of the response functions for inelastic electron scattering have been done 
by several authors ]13|] - p8[ . In the framework of Walecka's crw-model |l| the first calculations were performed for 
nuclear matter [|l3] - [|15[ and then applied to finite nuclei via the local density approximation (LDA) ||l3] , p^ . As a 
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result of many body correlations, in RPA the longitudinal response function is reduced by about 10-20% compared 
to the Mean Field (MFA) and Hartree approximation (HA). This reproduces at least qualitatively the trend of the 
experimental data. Similar results were obtained in a finite nucleus calculation ||l8[ . 

Despite the final curves of ||l3[ - [18| seem to confirm the same qualitative picture, they are obtained from two 
completely different approximation schemes. While vacuum polarization was neglected in ||l3] and pq] , it is included 
in the other calculations. In view of the results for the response functions one could argue that these seem not to be 
very sensitive with respect to Dirac sea effects. But this is not the case. The true situation is veiled by the fact that 
all calculations were performed only for the space-like kinematic region (q^ = Qq — Q^ < 0). The reason therefore is 
of course that only this regime is accessible by electron scattering. 

The aim of the present paper now is to draw attention to the fact, that effects in the space- and time-like energy 
regimes, which seem to be so well separated in view of the completely different experimental methods of investigation 
sketched above, nevertheless have an influence on each other. Especially it is our goal to show, that calculations which 
are performed to explain experimental data in one of the two kinematic regions have at least to be checked, whether 
for the other regime these provide results, which are physically consistent at least within the model used. So it will 
be shown that in the croj-model only calculations including Dirac sea effects will preserve causality. 

However, this becomes obvious only if one looks at the time-like part of the electro-nuclear response function. 
In RPA one would expect to find there a sharp peak corresponding to the a- and the w-meson respectively, which 
are the transmitter of the interaction in the Walecka model. Through their interaction with the nucleons these 
acquire the character of 'dressed' quasi-particles. This takes place in form of repeated creation and decay of virtual 
particle-hole and (if the Dirac sea is also taken into account) particle-antiparticle pairs which leads to a collective 
excitation of the whole nuclear matter system. Because a photon, emitted by an electron in a scattering process, 
interact with the nucleonic system via the same mechanism, these collective excitations manifest themselves as peaks 
in the corresponding response functions. More specific, there are four degrees of freedom: one scalar from the a and 
a longitudinal as well as two transverse from the w, where the latter are degenerated because of rotational invariance. 
Accordingly, two peaks are expected in the longitudinal response function (because of a mixing between the scalar 
and the longitudinal degree of freedom) and one in the transverse branch. These collective excitations in the time-like 
energy region have a direct back-coupling to the quasielastic bump in the space-like regime: they draw off strength, 
which causes the RPA-effect observed in former calculations [Q - |1^. In addition, these collective excitations 
should influence also the nucleonic properties which are studied in the heavy ion collision experiments mentioned 
at the beginning of this section. For example the masses of the mesons as well as the nucleons are shifted by the 



interaction |19[ which has consequences also for the subthreshold antiproton production |20|. This shows again how 
deeply connected phenomena in the two different kinematic regions are. 

The mechanism for the reduction of the quasielastic bump, described so far, of course makes sense only if the 
time-like collective excitations which draw off the strength from the space-like regime are physical degrees of freedom. 

As will be shown, this is only the case if Dirac sea effects are properly taken into account. A calculation where 
vacuum polarization is neglected shows only one meson peak in the longitudinal and no peak at all in the transverse 
branch of the electro-nuclear response function. Instead we find two peaks in the upper complex energy plane which 
means that these are not corresponding to well defined quasi particles but announce an instability of the nuclear 
matter system. Therefore vacuum polarization is crucial for the stability of nuclear matter. 

The purpose of the present paper is to study general properties of collective excitations in the nuclear matter system 
which are rather a consequence of many body effects then of the basic NN-potential actually used. Because it is not 
our intent to compare our curves with experimental data we use with the Walecka model a simple but nevertheless 
nontrivial NN-interaction. This has the advantage that for space- and time-like (real) energies the calculations may 
be done analytically which makes it easier to perform the numerical extension to complex energies. 

In section II the criteria for the stability of collective excitations will be discussed in general. In order to explore 
the upper complex energy plane with respect to destabilizing unphysical modes we have to perform an analytic 
continuation of the electro-nuclear response functions. This will be done in section III. Using dispersion relations we 
introduce a renormalization concept ||21[| which is new in this context. Instead of using dimensional regularization 
- as it was done in Mm and p7[ - we perform subtractions in the dispersion integrals for the polarization tensor in 
order to get rid of the infinities due to vacuum polarization. We show explicitly that both renormalization technics 
lead to the same results for the response functions which is an interesting result in itself. 

In section IV the response functions are at first discussed for real energies. In section V we look at the single 
particle properties in the space-like (quasielastic bump) and the time-like regime (particle-antiparticle excitations). 
In both cases as a consequence of RPA-correlations we find a reduction of the response functions compared to the 
independent particle model. This effect is most striking for the particle-antiparticle production amplitude which for 
high energies is reduced exactly by a factor 2. We show that this halving of the amplitude is a general result which 
holds true for any isospin independent interaction. This is one of the main results of the present paper. 

In section VI we discuss the role of collective excitations extending our considerations from real energies to the 



upper complex energy plane. Neglecting vacuum polarization we find the result mentioned already above, that the 
nuclear matter system becomes unstable because of excitations with complex energy z where Im(z) > 0. Finally our 
conclusions are given in section VII. 

II. STABILITY OF COLLECTIVE EXCITATIONS 

A possible mass change of light vector mesons inside nuclear matter as a consequence of medium modifications 
of the quark condensate < 9<z > (tI has recently drawn attention to calculations of timelike collective excitations in 
various models. An comprehensive overview of this work can be found in Q and the references given therein. Hatsuda 
et. al. estimate the mass reductions of p, w and <i> using QCD sum rules and compare these results with a calculation 
in the framework of the crw-model (which is suitable extended in order to include strangeness and the ^-meson) . In 
the latter approach which is similar to ours in this paper they obtain the effectives masses from zeroes of the inverse 
meson propagator in the long wavelength limes jfjl — > inside nuclear matter. 

Apart from this actual interest in density dependent changes of the quark condensate collective modes in nuclear 
matter were studied already previously by several groups. The first calculation within the crw-model was performed 
by Chin Ig] . Neglecting vacuum polarization he did an additional severe approximation [J evaluating the momentum 
integrals in the expression for the dielectric function. So Chin's results are partly in contradiction to the work 
of Lim and Horowitz 122] who calculated the density dependent momentum integrals exactly but neglected vacuum 
polarization too. These authors also addressed the problem of the stability of collective modes, however from a different 
point of view then ours in this paper. Lim and Horowitz investigated the density dependence of real collective zero 
energy modes (go = 0) in the space-like region. If for a given density (characterized by the Fermi momentum kp) and 
momentum \q\ such a mode arises, this indicates that in the uniform nuclear matter system an instability occurs which 
drives a phase transition to a spatial inhomogeneous state. A similar type of instability was studied by Furnstahl 
and Horowitz p3[ in connection with the so called Landau ghost which appears in meson propagators when vacuum 
polarization is included. Neither of these singularities will affect our calculations. The first type is excluded because 
we are only interested in nuclear matter at saturation density (kp = 1.4 /m~^) where the system is stable against 



zero energy modes (see fig. 12 of |22 ). The Landau ghost appears only for momentum transfers of \q\ > 2.55 m (m — 
nucleon mass) |23] which is much higher then the momenta we are considering n. 

In the following we examine the nuclear matter system with respect to a different kind of instabilities. In the previous 
calculations of collective nuclear matter excitations mentioned above, only the denominator of the meson propagator 
was analyzed on the real energy axes. Zeroes of this denominator are good candidates for corresponding poles in the 
propagator i.e. an excitation of the system. But, while even a zero of real and imaginary parts of the denominator 
is not always a sufficient condition for a pole (if the nominator also becomes small) , often only the vanishing of the 
real part for an given energy loq is already taken for an excitation. If the imaginary part of the denominator keeps 
sufficient small (and if the nominator is finite) this is correct and the propagator in the neighbourhood of cuq has 
the shape of a Lorentz curve corresponding to an excitation with a decay width depending on the magnitude of the 
imaginary part. 

In a paper by Jean et. al. ||2^ the denominator of the w-meson propagators has been analyzed similar to [g| and 
PI in order to find a mass shift as result of medium modifications in the aw-model. As it is one of the aims of the 
present paper, Jean et. al. also studied the influence of the Dirac sea on the effective lo mass comparing the full 
calculation with one where vacuum polarization is neglected and only Pauli blocking is taken into account. While 
in the presence of vacuum polarization all degrees of freedom are correctly reflected by a corresponding pole in the 
propagator they found in the latter case zeroes of the real part of the denominator simultaneously with a nonvanishing 
negative imaginary contribution. They concluded that "this nonzero imaginary part contributes to an unphysical (i.e. 
negative) decay width for the uj meson" . 

As mentioned above, at this point it becomes obvious that it is not sufficient to look only at the denominator of 
the propagator. In section VI we will show that the zero of the real part of the latter, found by Jean et. al., does 
not correspond to a collective excitation with "negative decay width" (whatever this means) at least not on the real 
energy axes. However, in order to find out what is really hidden behind this unphysical mode, it is necessary to 



'^As a result of this approximation, beside vacuum polarization also the density dependent Pauli blocking of NN excitations 
due to the filled Fermi sea is excluded from the calculation. 

^In a paper by Tanaka et al. |24|] it was shown that the Landau ghost is not a serious problem of relativistic many body theory 
because it can be eliminated in an appropriate way 



explore the full propagator (not only it's denominator) in the whole complex energy plane. As noticed already in the 
introduction, the neglect of vacuum polarization entails an instability of the nuclear matter system indicated by poles 
of the electro-nuclear response functions in the upper complex energy plane. 

In the case of inclusive (e,e')-scattering there are two independent response functions, a longitudinal and a trans- 
verse, which can be expressed by the corresponding matrix elements of the polarization tensor n^"^ {qo,q). If we choose 
a co-ordinate system were qis parallel to the a;-axes, Ii^'^{qo, q) can be written in the form: 



^^"{q^.q) 



I n^{qo,q) fnL(go,g) N^ 

fnL(go,<7) (f)'nL(qo,'7) 
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The Fourier transform of 11^*^(^0 7 ^'■ 
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is defined as the change 5j'^{x) of the nuclear current in x due to a change 5Au{y) of the vector potential in y, induced 
by the scattered electron: 

5f{x)^W''{x,y)5A,{y). (3) 

According to linear response theory pGl Il'^'^{x^y) can be expressed by the ground state expectation value of the 
commutator [j^{x)^ j'^ {y)]. Causality requires that H^'^[x,y) can be non zero only if x is in the forward light cone 
with respect to y, i.e. t^ > ty. Including a factor (-i) which is conventional, n^''(a;,y) can be written explicitly: 

Il^^x,y)^-i <0\[rix),f{y)]\0> 0(1, -ty). (4) 

Here |0 > denotes the ground state of the nuclear matter system consisting of a filled Fermi sea of nucleons above the 
Dirac vacuum of antiparticles. 

Now, for ty > tx the go^integral in (^ has to be closed in the upper complex energy plane. Then a zero contribution 
can only be obtained iiIl^'^{qo, q) is analytic there. On the contrary, if n^'^((7o, q) has poles in the upper half plane, this 
means that there is a response of the system preceding an external disturbance which is equivalent to an instability. 
Because in H , p2| , |23] only poles on the real go~axis were examined, non causal singularities could not be found by 
these authors. 

III. ANALYTIC CONTINUATION TO COMPLEX ENERGIES 

For nuclear matter (point-like Dirac particles) with an interaction mediated by a scalar and a vector meson (ctw- 
model m) the longitudinal and transverse electro-nuclear response functions within RPA are given by the following 
expressions Q 0: 

nW(9o,9')) = -nlAqo,q^) 



- e' 



(1 - J^Gv^vv)'~'s \^^sv) + 2Y^^'"^^^vv^sv + f5-(l ~ Gs^ss)Gy (2!!^ 



(1 — ^^t)n^^)(l — Ggllss) — ^G^Ggllj,^ 



(5) 



^Note that in contrast to (0) in ||l3] Il^^{x,y) is defined as the time ordered product 
— i < I T j'^{x),j'' {y) I >. For go > the Fourier transform of this expression is identical to our Tl^'^{qo,q} and differs only 
in the sign of the imaginary part for go < 0. 






(6) 



with: 



G.,„(g2) = , ^'," ^ . . (7) 

e,gs and g^ are the electromagnetic, scalar and vector meson coupling constant respectively. Tlab{(io,(i) {oL,b = s,v) 
is the Amplitude for a scalar or vector boson "a" to go into a nucleon-hole or nucleon-antinucleon pair which 
recombines to a boson "b" n. In former nuclear matter calculations m3 - |l3] Ilabicio,<l) was evaluated explicitly 
for real go employing Feynman rules. In order to get a finite vacuum contribution dimensional regularization was 
performed. In the following we use a method developed by Berestetskii in the context of QED pl| which is more 
appropriate for an analytic continuation of (^ and (|g) to complex energies z := qo +iqQ- We start from the imaginary 
parts of Ilaa for energy transfers qo > 0- Unitarity requires that these are given by the probability of boson "a" going 
into a nucleon-hole or a nucleon-antinucleon pair: 

Imn«,(go,9") = (Am*^ - q')Ioiqo,q'') +lmIlTi<l^rn*-') ; (8) 

ImUUqo,f) - -flo{qo,f) + hiqo, q') + g^Imn^ (-z', m*') ; (9) 

ImUUqo,q^) - -(4m*2 + q^)Io{qo,f) - ^hiqo,^) ~ <S,\vc,WZ{q\m*^) ■ (10) 

The density dependent contributions to these amplitudes are given by the integrals /„ over Fermi sea momenta up to 
the Fermi edge fcp: 

■^ k k-\-q 

+ (2ej-gor%0 -[££+,-+££])} . (11) 



Here ng — 8(fcF — |fc|) denotes the Fermi distribution of the mean field ground state and eg = vfc^ +m*'^ is the 

energy of a nucleon with wave vector k and effective mass m* . The (5-functions in the integrand of (O) make sure 
that the two terms give a contribution only in the space-like or the time-like region respectively. The momentum 
integral can be easily performed analytically [|l3|. 

The Dirac sea contributions of boson "a" going into a nucleon-antinucleon pair are given by: 



Im^:r(g^m") = "f^ (l - ^) ' e(g^ - 4m-) (12) 

and 

1 

1 / 4777*2 \ 2 / 0™*2\ 



Imnr(9^m*^) = -— (^1-— ^j (^l + __je((z^~4m*^) . (13) 

In addition to the "proper" transitions a — > (ph or NN) — > a {a — s,v) there is also a mixing of scalar and vector 
degrees of freedom. However, this is a purely density dependent effect. Because the vector current of the vacuum 
is zero, its commutator with the scalar current (which according to (||) enters the expression for the components 
Hst, of the polarization tensor) vanishes. The vector current of the nuclear matter system on the other side has a 



*For vector bosons ni,„ = IIJJ^ is a 4 x 4-tensor, the longitudinal and transverse matrix elements of which are denoted by nj;„ 
and n!J^. Because we have drawn out the coupling constants and isospin factors from the vertices, the analytical form of Tlw 
is the same for a;-mesons as for photons. A difference between Huiui, Hi^-y and Il-^-y will occur only in the numerical expressions 
of the vacuum contributions as a consequence of the renormalization procedure (see (pOf)). 



non-vanishing time-like component which is given by the energy density of the filled Fermi sea. The corresponding 
matrix element of the polarization tensor for qo > has the form: 

lmn,,{qo,f) = 2m*^hiqo,f) . (14) 

For qo < lmIlabiqo,q'^) is given by the negative of the expressions (||), (^), ( JU] ) and (14) respectively: 

lrRnab{qo,q^) = -lmUab{\qo\,q^) {qo < 0) , (15) 

i.e. the spectral functions linIlabiqo,q^) are antisymmetric with respect to go- 
Making use of the spectral representation: 



Tlab{z,q^) 



dqo lmIlabiqQ,q'^) 



qo 



(16) 



from the analytic expressions for the imaginary parts (H)-(tO|) and (p4) the full Ilab can be obtained in the whole 
complex energy plane. For go > 0, Imllab can be treated as a function of q^. Using the antisymmetry property (|l5| ) 
we can write dlq) also in the form: 



Ilab{z^,f) 



dqllmllab{qi,q^) 



q^-z^ 



(17) 



However (O) is directly applicable only to the density dependent portion of Hab because this has the appropriate 
asymptotic behaviour for go ~* oo- In the high energy regime the imaginary parts of Ilss and Ilw are given by the 
vacuum contributions Imn™'^((7^) and Imn"2^(q^), the asymptotic behaviour of which for go ~* oo,(q'^ ~ const.) 
follows from (|l|) and (^: 



Imn™^(g2,g-2)(xgg 
Imn™'=(go,g^) oc const. 



(18) 



Therefore in order to get from (UV) a finite result for the vacuum contributions, we have to perform in the dispersion 
integrals one subtraction with respect to gg for II^^, and two subtractions for H^j : 



uiTiz', f) = n:r (4„, f) + (z^ - 4 J 



m 



{z\t) = Knz'r;s,r) + {z'-zt,,) 



dql^Ton^^i^^q^ 

n (gg-z2.J(g2-z2) 



dli--%z^q^) 



Hz' - zlS^ 



dql 



dz^ 
lmli^,^%qlq^) 



T' (90-4^^(90-^^) 



(19) 



Using ([1^) it can be easily checked that the integrals in (|T^) are indeed finite. This procedure is equivalent to the 
corresponding subtraction of infinities in the dimensional regularization method p4[ , |l5[ or any other renormalization 
scheme. Accordingly, the finite values of n™"^, H™^ and dW!.'^'^ /dz^ have to be fixed at appropriate renormalization 
points z^. We have chosen to follow the "on mass shell" renormalization of llj] M: 



^There is a certain arbitrariness in the renormalization procedure for effective field theories because it can be argued that 
the masses of the scalar and vector bosons, mediating the interaction in these theories, do not have to be identical with the 
masses of physical mesons. Thus some authors |l3] , psi prefer to take e.g. q^ — Q as renormalization point . Fortunately the 
numerical results for the response functions are only slightly changed by this different choice. 



n™^(4^ = f + ml,m*=m)=Q 






dWZ{.z\m*=m) 



dn™=(z2,m* =m) 



dz^ 



= Sml 

-q^-\-m'^ 

= . (20) 



-q ^H-m^ 



Here Sm^ denotes a mass shift term which is introduced in order to keep the pole of the renormahzed w-propagator in 
the vacuum at the known mass of m^ = 783 MeV. With the boundary conditions ( |20| ) the (finite) vacuum contributions 
to the response functions ([19) are completely determined in the whole complex energy plane. 




IV. RESPONSE FUNCTIONS 

Having obtained explicit expressions for Habiz'^, 9^), the analytic continuation of the RPA-expressions 
can be easily obtained by the substitution qq —>■ z = qo + iq'^. On the contrary, the original expressions 
can be recovered from Hrpa(z, 9^) n by the prescription: 

nRPA(go, q^) = Hrpa(2: = go + «e, q^)- (21) 

Here it should be noted that the extrapolation to complex energies (O) is even necessary to obtain the correct result 
for the response function on the real energy axis. The reason is, that the ^-like peaks in ImHRPA which are caused 
by the collective excitations may be calculated only taking the proper limes according to (19|) (see section VI). 
The numerical values of the parameters of the theory are given in ^ : 

rris = 550 MeV m„ = 783 MeV 

92.87 (MFA) ^ _ ( 135.7 (MFA) *^^^) 

69.87 (HA) fft- - I 79, g (hA) 

The coupling constants are adjusted to the saturation property of nuclear matter at an energy of —16 MEV per 
nucleon for a density corresponding to a Fermi momentum of fcp — 280 MEV/c. While in mean field approximation 
(MFA) vacuum effects are neglected, the Hartree approximation includes these properly. The corresponding values 
for the effective mass of the nucleon are given by m* = 0.56 r??, (MFA) and m* — 0.78m (HA) respectively. 

Performing the limiting procedure according to {Uw the electro-nuclear response functions can be obtained from 
H(z,g2): 

Rl,t('?o, g") = IninL.T(go + «e, f)- (23) 

TT 



V. SINGLE PARTICLE PROPERTIES FOR SPACE- AND TIMELIKE MOMENTUM TRANSFERS 

Part a) of figs. 1-4 show Rl and Rt for a momentum transfer of \q\ = 1.5 fcp which was a typical value in former 
electron scattering experiments ]1^ , ||ll| . Part b) of the figures show the corresponding real portions Re Hl and Re Ht 
of the longitudinal and transverse matrix elements of the polarization tensor (see equ. (0)). The energy ranges from 
the space-like region (go < \q\) up to the deep time-like regime (go > \<l\)- While vacuum polarization is included in 
figs. 3 and 4, in figs. 1 and 2 purely density dependent contributions arc shown. As mentioned already in section I, our 
goal is to point out qualitative differences between these two approximations rather then to compare the results with 
experimental data (which for the quasielastic regime - where such data are only available up to now - has been done 
by several authors [1^, |l^, Q, Q). Therefore we did not include form factors but performed the calculations for 



^The same holds true of course for Ilab according to (|l9|) 



point-like Fermions in nuclear matter. The dashed curves show the results including only one polarization bubble 
(corresponding to the first term in equs. (|5|) and (|6|) respectively), while in the full curves RPA-correlations are taken 

into account according to the the second expressions in (0) and (p). The structure on the left sides (w/fcp ~ 1) of the 
four figures represents the quasielastic bump which is measured in (e,e') scattering experiments. According to ( |2^ ) 
the effective mass m* is higher when vacuum contributions are included. This causes the quasielastic peak to be a 
bit narrower in this case [Q, fl^l. The cusp at w/fcp^O.lS in the RPA-curve of fig. la is due to an almost zero point 
of the denominator in equ. (g|) which indicates that long-range correlations are important. 

From figs, la and 2a one easily reads off that the purely density dependent contributions to the response functions 
are zero for lu -^ oo. Therefore RenL,T in figs, lb and 2b can be obtained from the curves in Figs, la and 2a directly 
(i.e. without any subtractions) by means of equs. ( [l7| ) and (p3|). The quasielastic bump in the response functions 
has its counterpart in a 'oscillator-structure' in the real parts of IIl.t- The cusp at w/fcp^O.lS in the RPA-curve of 
fig. la corresponds to a small 'shoulder' at the same energy in fig. lb. 

With vacuum polarization as well as without, many body effects which are taken into consideration in the RPA 
provide a considerable reduction of the longitudinal and a smaller decrease of the transverse response functions. The 
reason for this suppression in the space-like region has its origin in the time-like regime and consequently keeps 
hidden for electron scattering experiments: strength is transfered to collective modes which in the RPA-curves appear 
as sharp peaks in the so called 'unphysical region' \q\ < uj < y/4ni*'^ + q"^. These peaks will be discussed in more 
detail below. Despite the RPA-curves are quite different in shape, the omission of vacuum polarization does not 
change drastically the qualitative picture in the quasi-free regime. Comparing figs. 1 and 2 with figs. 3 and 4 one sees 
that the last statement holds true for the response functions itself as well as for the the real parts of IIl x- 

This changes completely in the time-like energy range. For qq > ^/4m*^ + q^ the production of NN-pairs from the 
Dirac sea becomes possible. According to ^ and ( [l^ ) for the longitudinal response (fig. 3a) this vacuum polarization 
has a "square root-like" go-dependence (but the slope keeps finite) at threshold and approaches a constant value for 
increasing energies, [j 

From equs. (Bh and (nfl) one finds for qo >> \q\ as a relation between the longitudinal and transverse matrix 
elements of Im H^^ : 

ImnT (q2^^2^ g2 g2 



ImnL^((j2,(f2) g-2 ^ 



■2 



As will be shown below, for high energies not only the one bubble contribution but also the RPA-corrections to 11 ' 



in (^) and (^) are proportional to Hy^^ and consequently also obey the relation (g^ Q Hence vacuum contribution to 
Rt in fig. 4a in comparison with Rl in fig. 3a is weighted by an additional factor g^ which causes a rapid increase for 
energies above threshold. So in order to show the low and high energy branches in the same figure we had to scale 
the transverse NN-contribution by a factor 0.05 . 

In the presence of nuclear matter the production of particle-hole pairs from the Dirac sea is inhibited if the final 
state has momentum lower then kp. So the Fermi sea causes a reduction of the vacuum amplitude for the process 
e+e" -^ NN. Figs, la and 2a show this negative density dependent contribution to the longitudinal and transverse 
response function in the time-like energy regime. The big slope of Rl,t at threshold is reflected as a cusp structure 
in the real parts of IIl^t (see figs, lb and 2b). Again the transverse branch is weighted by an additional factor 
g2 (compared to the longitudinal contribution: see equs. (|9|) and (|lO[)) which results in an enhancement of the 
corresponding time-like counterpart of the quasielastic peak. 

From figs, la and 2a it is easily read off that approximations which neglect vacuum polarization suffer from the 
deficiency that the response function is no more positive definite over the whole energy range. The sum of vacuum 
and density dependent contributions (figs. 3a and 4a) however is always positive ^. As noticed already in sec. Ill, with 



'^Note that the curves in figs. 3a and 4a show the sum of vacuum polarization and the corresponding density dependent 
contribution. But it follows from (H) - ( |ll| ) (and can be explicitly read oflt from figs, la and 2a) that in the time-like region the 
latter are nonzero only for the small energy range min(2e|q-|/2, eg ^i^i + efcp) ^ 1o < ^E +|g-| + ^fcp- So for higher energies in the 
presence of nuclear matter the shape of the curves is the same as for the vacuum. 

*From the solid and dashed curves in figs. 3 and 4 one sees already that in the high energy limes the two approximations can 
differ only by a multiplicative factor which - as will be proven analytically below - has the value 1/2. Therefore in the following 
considerations on the asymptotic behaviour of the polarization propagator we do not have to distinguish between IIjj^p^ and 

®If one compares figs, la and 2a with figs. 3a and 4a it seems that the (negative) time-fike density dependent contributions 



regard to vacuum polarization IV^^ does not fall off to zero for qq -^ oo but approaches a constant value (see equ. 
(|l8|)). As a consequence of our remarks above this holds true also for the longitudinal response function (see fig. 3a) 



Therefore the real part of IIl in fig. 3b cannot be obtained directly from the curve in fig. 3a making use of (17) as 
it was the case in fig 1. Instead, according to the first equation of (|lj) one subtraction has to be performed. The 
meaning of this renormalization procedure becomes more obvious, if one splits up the vacuum contribution in fig. 3a 
in a constant part (which starts at the NN-threshold) and a contribution which has an asymptotic behaviour like 
1/(7q for qo ^ oo (compare equ. (p^). More specifically we write RL(qoi'?*^) i^ the form: 

RLiqlf) = RL(go' - ^, f) e((z' - 4m* 2) + ARL(go', -f") (25) 

where by construction AKi,{qQ,q^) provides a finite contribution to the integral ( p7| ) which is harmless. The lead- 
ing order high energy contribution to Re IIl comes from the 0~function which may be integrated performing one 
subtraction: 

RenL(gg,a--9o^ / —^JT^ -.-HqI)- (26) 

J "' 9o wo %) 

4m' '^+q^ 



The negative sign of the principal value integral in (|2g) follows from the definition ( |23| ) of the response function: 
Rl ~ — ImllL > 0. In this way we obtain the same result as by using dimensional regularization uM, |17[ |. 

Because of ( |24| ) two subtractions are needed in order to calculate RenT(goi Q^) from the response function in fig. 4a 
(i.e. one has to integrate KT{qQ,q'^)/qo and multiply the result with q^). Analogous to ( p6| ) one obtains to leading 
order in gp: 

RenT(go','?")-^'?oM'?o)- (27) 

This result can be easily checked in fig. 4b. Hence the relation (p4) is valid also for the real parts of nL,T, i-e. in the 
high energy limit we get the final relation: 

nT{qlf)^q^Ill,iqlq^) (28) 

which will be used in the following F^. 

RPA-correlations lead to a dramatical decrease of NN pair production by a factor 1/2. In our approximation 
of point-like nucleons without anomalous magnetic moment this suppression is a consequence of isospin symmetry. 
For the transverse response function this can be read off most easily from equ. (|6|). Because the photon in our 
approximation couples only to protons while meson exchange is charge-independent, the 'electromagnetic bubbles' 
H'^y for single particle-hole pair creation as well as in the numerator of the RPA-correction are weighted by a factor 1/2 
in contrast to the 'mesonic' bubble in the denominator. For high energies the g^-dependent part of this denominator 
becomes much larger than 1. This can be checked substituting 11,^^ according to (P8|): the leading factor q^ cancels 
the 1/g^-decrease of the propagator Gy{q'^), so that the RPA-dcnominator in (pi) goes like {l~ g^ ^vv) ^ ^~ st ^^{lo) 
(see equ. (|2q)). Hence for high energies the 1 in the RPA-denominator may be neglected. Then the rest of the 
denominator in (0) cancels against the corresponding part in the nominator so that (because of the isospin-factor 
1/2) one gets from the RPA-correction a contribution which is exactly (-1/2) times the single NN-bubble contribution 
PI. The same holds true for the longitudinal part (S) of the polarization propagator. As pointed out in the section 



preceding equ. (14), H^^, is a purely density dependent contribution. From equ. (ll|) it follows (and can be easily 
checked in figs. 1 and 2) that Hst,(go,9^) = for q^ > er , i-i -f ea-f- Consequently only the third term in the nominator 
and the first term in the denominator of (0) are nonzero for go high enough. So in this limes Hj^p^ can be written: 



start at qo ~ 4feF which is about 1 fcp below the threshold for pair production. But one has to keep in mind, that these 
curves show the result of two different approximations. As was pointed out in section IV the effective nucleon mass is higher 
when vacuum polarization is taken into account properly. Consequently the NN-threshold min(2e|q-|/2, ^% _\^ -\- ^k-p) is shifted 
according to the corresponding mass shift from m* = 0.56 m to tti* = 0.78 tti. So in figs. 3a and 4a the density dependent 
contributions start at the same threshold as vacuum polarization 

^ Here of course we recover the tensor structure: llJJac = [s'"' — Q^g"/? ]Rvac{q )■ So equ. (|28|) proves that our renormalization 
procedure preserves current conservation. 

'^^Note that the starting energy of this asymptotic behaviour depends on the numerical value of gl which is a parameter of the 
model. If gl » 1 (compare equ. (P2)) the reduction by a factor 1/2 starts practically already at threshold which is confirmed 
by figs. 3 and 4. 



[nRPAlL (*, f)) - y n^„(go, f) + e^ J^;!^;^, . (29) 



Completely analogous to the RPA-correction to IIt , the 1 in the denominator of (|29|) may be neglected and one gets 
a negative contribution of just one half of the single pair production amplitude. 

This argumentation on the one hand is very general because it is based only on the isospin symmetry of the NN- 
interaction and holds true beyond the crw-model. So for point -like nucleons we would expect a suppression of the 
cross section e'^e~ —^ pp by 50%, inside nuclear matter as well as in the vacuum. On the other hand we know that 
nucleons are not point-like particles, which may invalidate (at least partly) our reasoning. Nevertheless, we think that 
our model calculation has shown that many body correlations can produce effects far beyond the usually assumed 
level of 10-20%. This result should not alter dramatically taking the finite size of nucleons properly into account. 

An observable which would be directly affected by such many body effects is the ratio: 

R = "(f f ^ ^"^'^""^^ (30) 

Within the simple quark-parton model (QPM) - which corresponds to the one bubble approximation in our calcula- 
tions ~ R is given by: 

R = 3^e2, (31) 

where the index q runs over all quark flavours with charge Cq. A review on the value of R measured in several 
experiments can be found in ||27| . Unfortunately our considerations concerning the process e^e~ — > pp cover only a 
fraction of R which makes it difficult to compare directly with these experiments. Nevertheless, a thorough analysis 
of the above general isospin argument on the quark level seems to be worthwhile but goes beyond the scope of the 
present paper. 

VI. COLLECTIVE PHENOMENA 

Another remarkable difference between the full and dashed curves in figs. 1-4 are the narrow peaks in the time-like 
energy regime below the threshold for pair production. These belong to the collective excitations already discussed 
in sections I and II. The correct structure of these poles in the real and imaginary part of the polarization propagator 
may be obtained only if the calculation is extended to complex energies. The sharp 5-peaks in the response functions 
(see part a of figs. 1-4) are the result of an analytic continuation back to the real energy axis according to (P3). 
For Im(2;) = the 'one bubble' contribution to the response function is zero in the unphysical region \q\ < qo ^ 
min(2e|g^/2i £fc _m + ^fcr) 0- ^^ ^ consequence, also the nominator of the RPA-corrections in equs. (ph and (0) 
vanishes for these energies. Hence zeros of the denominator, which are the signal of collective excitations, lead to an 
indefinite expression for the RPA-formulas. 

Former calculations of the electronuclear response functions [Q- [Q were restricted only to real spacelike energies 
(more specific: to the quasi elastic bump) and therefore did not have to care about this problem. On the other hand, 
the authors who looked for timelike excitations 0, l|], [g2|, [g5[ only calculated the denominators of equs. (||) and (||) 
and looked for zeroes in this quantity. As will be shown below, for real excitations this method provides the correct 
spectrum, but in the case of an instability such a treatment may be misleading. 

From counting of degrees of freedom (see section I) one would expect two peaks in the longitudinal and one (twofold 
degenerate) peak in the transverse response function. But figs. 1 and 2, where only density dependent contributions are 
taken into account, show one single excitation in the longitudinal and none at all in the transverse branch. However, 
including vacuum polarization the curves in figs. 3 and 4 show the correct number of peaks. 

Here it is instructive to look at the RPA-dcnominators in both approximations which are shown in fig. 5 for the 
longitudinal branch (the following arguments hold true also for the transverse response). The large peaks in the real 
parts of both figures come from the free cr and ui propagator in the denominator of equ. (ra) with masses given in equ. 
(p2|). The dispersion of the dressed mesons inside nuclear matter is determined by the zeroes of the full curves in 



^ There are no single particle excitations in this energy range 
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fig. 5 (but, as said before, only if the imaginary part, i.e. the dashed curve in fig. 5, vanishes too). Including vacuum 
polarization (fig. 5b) we find a reduction of the effective meson masses because the zeroes are situated below the 
peaks of the free propagators. This result is in accordance with previous calculations |g[, Q, |25|. On the contrary. 



neglecting Dirac sea effects, the real zeroes of the RPA-denominator (from which only the first corresponds to a well 
defined excitation) are situated above the free meson poles which means an enhanced mass of the quasiparticlcs. 

In order to understand the reason for this difference in the two approximations it is helpful to take a closer look at 
the shape of the peaks in fig. 5. If one compares fig. 5a with fig. 5b it is easy to verify a different sign of the meson 
poles. While in fig. 5a the real part of the denominator is positive below the poles and negative above, the zeroes 
corresponding to the quasiparticlcs also lie above the poles. Exactly opposite is the situation for the calculation 
including vacuum polarization (fig. 5b) which leads to the reduction of the meson mass. 

Now, the mechanism how Dirac sea effects cause this negative mass shift becomes clearer, if one looks at the 
denominator of the RPA correction in equ. (|^). In the energy range between the photon point q^ — and the 
threshold for pair production q'^ — m*^ the mixing term of scalar and vector degrees of freedom 11^^ can be neglected 
in comparison with the pure terms H^^ and Hgs- So for these energies the RPA denominator may be approximated 
by the expression (1 — ^GyII\^^){l — Gs^ss)- The poles of the propagators Gs and Gy provide the peaks discussed 

above, the signs of which depend on the sign of Hy^ and Uss in the neighbourhood of these poles respectively. Because 
without RPA corrections IlJ^^ is directly proportional to the longitudinal response function (see the first term in equ. 
(g)) its behaviour can be red off from the dashed curves in figs. 1 and 3. For the energies under consideration the 
imaginary parts vanish. Neglecting the Dirac sea, the dashed curve in fig. lb for q"^ < q^ < q'^ + m*'^ has positive sign. 
Because the scalar and vector meson masses are well separated, near the vector pole (1 — G^IIss) ~ 1 > 0. This leads 
to the peak structure described above and an enhancement of the to mass inside nuclear matter in this approximation. 

Comparing fig. lb with fig. 3b it is easy to see, that vacuum polarization provides a negative contribution to the 
real part of Hyy . The result of this is the switch of sign in the free cu pole observed in fig. 5b with the consequence 
of a reduced mass of the dressed quasiparticle. According to equ. (nS) the magnitude of the negative Dirac sea 
contribution and thus also the meson mass shift depends on the amplitude for nucleon-antinucleon pair production. 
The threshold for this process inside nuclear matter is determined by the effective nucleon mass m* which is the free 
nuclcon mass m reduced by a term proportional to the scalar density < -ip^ > [|l|. Hence analogous to the results 
of the QCD inspired approach M , with a meson mass shift due to medium modifications of the quark condensate 
< uu > also in the Walecaka model a scalar density is responsible for the changes in the dispersion of the mesonic 
degrees of freedom. 

Despite the preceding considerations, as an illustration of our method, gave us a qualitative picture of the underlying 
mechanism for the meson mass shift, it is not the purpose of the present paper to analyze this effect in more detail. 
Rather it is our goal to point out the importance of vacuum polarization for the stability of the nuclear matter system. 
In |g5[ Jean et. al. argue that, neglecting vacuum polarization, the zero in the real part of fig. 5a at w/fcF~4.5 belongs 
to a collective excitation with "an unphysical (i.e. negative) decay width" because of the nonzero negative imaginary 
part of the RPA-denominator in this energy range. However, according to figs. 1 and 2 for real energies the response 
functions itself do not show any sign of such a unphysical mode. The cusp like structures at ui/kp^A in the real parts 
of the longitudinal and transverse matrix elements of the polarization tensor (figs, lb and 2b) are already present in 
the one bubble approximation and therefore are not the result of collective phenomena. 

The true situation unfolds itself only if the calculation of the response functions is extended to complex energies (see 
figs. 6 and 7). The peaks which are missed on the real energy axes reappear as undamped poles in the upper complex 
half plane. So in contrast to the prediction of [ p5[ an approximation which neglects vacuum polarization does not 
lead to a negative decay width of the lu meson but to collective excitations which are unphysical because, according 
to the discussion in section II, they manifestly destroy causality. Consequently calculation within this approximation 
cannot provide reliable results. 

VII. CONCLUSIONS 

We studied qualitative properties of the electro-nuclear response functions for nuclear matter with respect to 
single particle and collective excitations. Performing the calculations exemplary in Walecka's aw-model, the main 
results were obtained from general arguments which should hold true far beyond this simple model for the medium 
NN-interaction. 

In contrast to former calculations where the main attention was concentrated on the quasielastic region, we extended 
our considerations 1) to the time-like energy regime and 2) to complex energies. The latter was achieved making use 
of dispersion relations. As an alternative to dimensional regularization we introduced a renormalization scheme which 
is more appropriate to our approach: infinite vacuum contributions were removed performing suitable subtractions 

11 



in the dispersion integrals. With the appropriate boundary conditions this method leads to the same results as 
dimensional regularization. 

Neglecting vacuum polarization we found that there is a violation of causality because of unstable longitudinal and 
transverse collective excitations in the upper complex energy plane. Including vacuum polarization, these collective 
modes are shifted to the real axes and change to proper excitations of the system. So we have shown that vacuum 
polarization is crucial for the stability of nuclear matter. 

In addition we have shown that RPA-corrcctions reduce the production amplitude oi pp pairs by a factor 1/2. This 
result is remarkable because it holds true for any interaction between point -like nucleons which preserves isospin 
symmetry. Our argument may have interesting consequences for the process e'^e~ —>■ pp. 
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FIG. 1. (a) Longitudinal response function and (b) real part of the longitudinal matrix element IIl of the polarization tensor 
for nuclear matter at a momentum transfer of 1.5 fcp without vacuum polarization. The dashed curve shows a calculation for 
space-like and time-like energies neglecting RPA-correlations, while these are included in the full curve. 
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FIG. 2. Same as fig. 1 for (a) the transverse response function and (b) the real part of the transverse matrix element IIt of 
the polarization tensor. 
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FIG. 3. Same as fig. 1 including vacuum polarization. 
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FIG. 4. Same as fig. 2 including vacuum polarization. Note tliat in (a) because of the fast increase of the transverse vacuum 
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FIG. 5. Real part (full curves) and imaginary part (dashed) of the longitudinal RPA-denominator according to equ. (@). 
While in (a) vacuum polarization is neglected in (b) these contributions are included properly. 



17 





/ 


" 7 




K' 


/ 
/ 




0.8 - 
0.6 - 
0.4 - 


'Wiljim^^M 






/ 




% 






1 




f 


WmSSm^M 




0.2 - 


MMwWmM ' 








/ 






- 


r ■^^ 


i I 

I 1 

5 : 














-0.2 - 

-0.4 - 


/ / 






-—- 


- 



















^^y 


-0.6 - 


/ / 
/ 










y 0.8 ^. 


1 ■ ' 2 " '^ 

Re(z)/ 


% 


4 


5 


6 


/ 0.4 v:y^ 

0.2 \<^ 



FIG. 6. Extension of the longitudinal RPA-response function to complex energies neglecting vacuum polarization (compare 
the full curve in fig. la). The big singularity in the upper complex energy plane is the signal for the violation of causality in 
this approximation. 
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FIG. 7. Same as fig. 6 for the transverse response function. 
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